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Abstract 



We use the method of discrete loop equations to calculate exact correlation functions 
. of spin and disorder operators on the sphere and on the boundary of a disk in the 

c = 1/2 string, both in the Ising and dual Ising matrix model formulations. For both 
the Ising and dual Ising theories the results on the sphere are in agreement with the 
KPZ/DDK scaling predictions based on Liouville theory; the results on the disk agree 
with the scaling predictions of Martinec, Moore, and Seiberg for boundary operators. 
Q^' The calculation of Ising disorder correlations on the sphere requires the use of boundary 

1^ ■ variables which have no matrix model analog. A subtlety in the calculation on the disk 

arises because the expansions of the correlation functions have leading singular terms 
which are nonuniversal; we show that this issue may be resolved by using separate 
^ ■ cosmological constants for each boundary domain. These results give evidence that 

the Kramers- Wannier duality symmetry of the c = 1/2 conformal field theory survives 
coupling to quantum gravity, implying a duality symmetry of the c = 1/2 string even 
in the presence of boundary operators. 
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1 Introduction 



Our understanding of 2D Euclidean quantum gravity coupled to conformal matter fields has 
improved dramatically in the last decade, largely due to the development of matrix model 
methods. The connection between the discretized and continuum formulations of the theory, 
however, still involves a number of unresolved issues. Although the relationship between bulk 
operators in the continuum and in matrix models has been fairly well established |^, |^, Q 
(see for a review), the boundary operators in these theories are less well understood 
and have a more subtle relationship |^, ^, ^. While the matrix model formalism is 
naturally suited to the calculation of correlation functions of boundary operators 0, few 
explicit computations of such correlation functions between nontrivial matter fields have 
been performed (one calculation of this type is described in 0). 

In this paper we consider the simplest theory with nontrivial matter fields, the c = 1/2 
model which describes a single free fermion coupled to Liouville gravity. We compute exact 
correlation functions of operators in the discretized theory, both on the sphere and with op- 
erators inserted on the boundary of a disk. The c = 1/2 theory has two distinct descriptions 
as a matrix model, which can be constructed by coupling Ising spins to a triangulation or its 
dual graph respectively |[Ty, |TT], |r^. Although Kramers- Wannier duality relates these two 
formulations of the Ising model in fiat space, it is not clear a priori that this duality still 
holds after coupling to quantum gravity. In fact, it has been pointed out that there is a 
puzzle if one assumes that this duality exists |13|, in that the leading singular term in the 
two-disorder-operator correlation function on the boundary of a disk in the standard c = 1/2 
matrix model seems to disagree with results for the scaling dimension of the disorder operator 
in the dual theory. We resolve this puzzle by calculating the disorder operator correlation 
functions in both formulations of the theory, and showing that the universal behavior of 
these correlation functions is identical in the two models. On the disk, the leading singular 
terms in the two models agree, but are nonuniversal; the universal terms are determined 
by using separate cosmological constants for each boundary domain. The results presented 
here give strong evidence that Kramers- Wannier duality persists under coupling to quantum 
gravity. 

One particularly interesting feature of the calculation on the sphere is that, although the 
scaling behavior of the spin operator in the usual matrix model formulation of the theory can 
be determined by calculating the correlation function of the disorder operator in the dual 
formulation |T^, it does not seem to be possible to calculate the correlation function of the 
disorder operator (the spin operator in the dual formulation) by matrix model methods. To 
calculate this correlation function we use the geometric loop equation methods detailed in 



14 1 , taking advantage of the freedom in that formalism to use different boundary variables 



than those associated with the matrix model description. 

In Section 2, we review the matrix model descriptions of the c = 1/2 string as an Ising 
model coupled to dynamically triangulated 2D gravity, both in the standard and dual for- 
mulations. We include a discussion of the cubic equations satisfied by the homogeneous 
disk amplitude in each of these models. In Section 3, we take the continuum limit of the 
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homogeneous disk amplitude in both models, and compare to the predictions of Liouville 
theory. In Section 4, we compute the correlation functions of spin and disorder operators on 
the sphere in the Ising and dual models. In Section 5, we derive the correlation functions on 
the disk and compare with Liouville predictions. Section 6 contains concluding remarks. 
While this manuscript was being completed, we received a preprint by Sugino and Yoneya 



15| containing some results which overlap with those presented in Section 5. 



2 Discrete c = 1/2 models 



There are two ways of describing the c = 1/2 string in terms of a matrix model, each of which 
has been analyzed extensively in the literature. The original formulation of this theory as 
a matrix model is implemented by putting an Ising spin at the center of each triangle in a 
dynamically triangulated gravity theory |ll| . The corresponding matrix model action is 
given by 



Z = JbUBV exp 



-Tr f/2 + -Tr 



cTr UV 



3 



(1^ 



where g and c are coupling constants and U and V are N x N hermitian matrices. (Note 
that the coupling constant c is unrelated to the central charge c = 1/2. An unfortunate 
historical coincidence of notation makes it difficult to change either symbol; in this paper 
the symbol c always refers to the coupling constant except when it appears in the equation 
c = l/2.) 

In the alternate formulation of the c = 1/2 theory, Ising spins are coupled to the vertices 
of the dynamical triangulation. This corresponds to the Kramers- Wannier dual to ([T|), and 
may be thought of as a special case of the 0(n) model on a random surface [|1^, |16|, |1^. The 
action for this matrix model is 



Z = JdXBY exp (^-N 



[I - C) n (1 + C) ^ 



9 



(Tr + 3Tr XV'^ 



(2) 



In this model, edges of the dual graph connecting two Y matrices (the Y propagators) 
correspond to boundaries of constant spin domains in the dual Ising theory. In the remainder 
of this paper, for simplicity, we refer to the two models (|^) and (H) as the "Ising theory" and 
the "dual theory" respectively. 



2.1 Duality 

The two matrix model formulations of the c = 1/2 string are ontologically quite different; 
in the Ising theory spins are associated with faces of the random triangulation, whereas in 
the dual theory spins are associated with vertices of the dual graph of such a triangulation, 
whose coordination number is unconstrained. It is therefore interesting to recognize that a 
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simple coordinate transformation W7 



X 
Y 
9 



9 



(3) 



reveals that the partition functions for the two models are identical in the absence of bound- 
aries. The equality of the partition functions for these two models, however, does not imply 
that the Kramers- Wannier duality of the fiat space theory survives the coupling to quantum 
gravity. The transformation has a nontrivial action both on the operators and the states 
of the theory, taking the spin operator U — V in the original theory to the disorder operator 
Y in the dual theory, and taking free boundary conditions in the original theory to fixed 
spin boundary conditions in the dual theory (and vice- versa). To demonstrate that duality 
persists after coupling to quantum gravity, on the other hand, it is necessary to show that 
the correlation functions of the theory are the same in both formulations. In fact, we will 
show in this paper that unlike the partition function, the correlation functions in the two 
theories are not identical at the discrete level, but that the universal parts of the correlation 
functions which correspond to the continuum limit of these models are identical. 



2.2 Solution of models 

Because the action for the Ising theory is of the form 

S{U, V) = Tr ifuiU) + MV) + UV) 



(4) 



the matrices U and V can be simultaneously diagonalized []T9|, |I0|, |Tl[ , and the method 
of orthogonal polynomials can be used to solve for the correlation functions 



Pn 



— (Trt/") . 



(5) 



These correlation functions describe the disk amplitude with a boundary condition corre- 
sponding to all Ising spins being fixed in the same direction. Using loop equation methods 
it can be shown that the generating function 



1 



(6) 



fc=0 



satisfies a cubic equation. (We evaluate all matrix model expectation values in the large- 
limit.) This cubic equation has been derived by a number of authors using different methods 



m, H, |13|, H, im and is given by 



/30' + /2</)' + /l</' + /o 



(7) 
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where 



fo 

h 
h 



-g^ + 9^(2- gpi) u - g {1 + c - 2gpi + ^^2) 

- (5(2 + 51(1 + c)pi - 2g'^p2 + g^P3 - c(l - c^)) 
g^ — 2g'^u + g{l + c)u'^ — (c(l — c^) + g'^) + g {I — gpi) 
u^{2g'^ -2gu + cu^) 
gu^ . 



The coefficients Pi,P2,P3 can be determined either from the analyticity properties of this 
cubic, or by the method of orthogonal polynomials, pi and ps are given as functions of c and 
g by the relations 



Pi 



(3^ [(1 - 2(3f (-1 + 4/3 - 3/32) + (1 - 8/3 + 17/32 - 12/33) c + c^-c^ 



4(l-2/3)^^?3 



(9) 



and 



P3 



8(1 - 2/3)Ml - 5/3 + 6/3' - 2/3^ 



16(l-2/3r^;5 

+2(1 - 2/3)^ (-9 + 77/3 - 209/3^ + 235/3^ - 90/3^) c 
+ (6 - 131/3 + 658/3^ - 1390/3^ + 1364/3^ - 531/3^) c' 
+2 (6 - 2/3 - 71/32 ^ ^35^3 _ gQ^4\ ^3 



(10) 



+6 -2 + 3/3 + 5/32 - g/?'^ + 6 (1 - /3) + (-2 + /3) 



where /3 satisfies the quintic equation 

(/3 + c - 1)' + (1 - 2/3) [/3 (/3 + c - 1) (c^ - (1 - 2/3)') - 2(7^ (1 _ 2/3) 
The relation 



1-c 

P2 = Pi 



(11) 

(12) 



gives P2 in terms of pi |T^ . 

Although the dual theory is directly related to the Ising theory through the coordinate 
transformation (^), the dual theory is somewhat more difficult to solve by matrix model 
methods since the action is not of the form (^). Nonetheless, because the action is quadratic 
in y, this matrix can be integrated out, allowing the model to be solved [|r^, [T^, [T^. 

In the dual theory, the correlation functions 



Pn^^ (Tr X") 



(13) 
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correspond to the disk amplitude with boundary spins which are all identical, but which 
may be of either spin. Note that these correlation functions also describe the Ising theory 
with free boundary conditions, since X corresponds to (U + V)/\/2. Inserting a Y into the 
correlation function corresponds to fixing the dual spins on the boundary to be of opposite 
values on the two vertices separated by the edge labeled by Y. Thus, each Y corresponds to 
an insertion of the disorder operator. 
The generating function 

_ 1 / 1 \ °° 



for the dual theory satisfies a cubic equation similar to (|3) '^J} This cubic equation is 

/30' + /20' + /l0 + /o = O , (15) 

with coefficients 



(16) 



/o = -4c^^(l + xpi + x'^p2) + 2cgx{3 - c)(l + pix) - 2cx^(l - c^) - g'^x'^ 

= 4cf - 2cgx{3 - c) + 2cx^{l - c^) + gx^{l - 3c) - 2fpix^ 
/2 = x^{g^-gx{l-5c)-2cx^{l + c)) 

h = gx^ ■ 

The coefficients pi,p2 can be related through (^ to the coefficients pi 

Pi = V2pi (17) 
^ _ {1 - c^)pi - g{l + gps) 

P2 — • 

3 Disk amplitudes 

We now take the cubic equations satisfied by the homogeneous disk amplitudes in the Ising 
and dual theories, and expand around the critical point to extract the continuum limits. 

3.1 Ising theory disk amplitude 

The critical values of g and c can be determined by finding the radii of convergence for f3. 
They are given by [|10], ^ 

^ ' -1 + 2^7) (18) 



27 
and 

g, = 3-^/^Vl0[-l + 2V7)~' . (19) 
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From Eq. (|TT|) it follows, using and (|I^), that 

(3c = li5-V7), (20) 

and therefore that 

9cPic = k^-V7) (21) 
5 

and 

^cP3c = Y^(-699 + 40-7i) . (22) 

We can now compute the critical expansion of /3, pi and p^. The expansion for g, written as 

g = 5(cexp(-e^t) , (23) 

is expressed in terms of e^, where the continuum limit corresponds to taking e — 0, with 
e scaling as an inverse length. This choice is natural since the parameter t is conjugate to 
area. From this expansion it follows that 



22 + IOV7' 



gpi = ^ — ^(1-^ '-eH 



36 216 



and 



-699 + 40-71 4(121 -5- ri 



= ""100" " ^ 25 ^ 

^9j|^^8/3^4/3 ^ 3j|i ^10/3^5/3 ^ ^^^4^2) _ (25) 

Turning now to the expansion of 0, we take 

u = Ucexp{—ez) = gd^A — V7^ exp(— ez) , (26) 

so that z is the boundary cosmological constant in the continuum limit, conjugate to bound- 
ary length. The critical value of u is derived from the condition that at u = the cubic 
has repeated roots. Taking the expansions (^), (0) and ( pSf ) for g, gpi and gp^ above and 
substituting into (0), we find [|13|, ^ 



3(2 + 77) / (l + 2v^' 



ez 



+ e^/=^0(4/3)(^,t)+e^/'0(5/3)(^,t) + O(e') , (27) 



10 \ 3 

where 0(4/3)(z, t) solves the equation 

270^4/3) - 9 ■ 5^0(4/3)^^^^ - SOt^ + 20(2 + Vrftz'^ - (2 + Vt^z^ = , (28) 

7 



and (^(5/3) is given by 



<P{5/3){z,t) 



20(4/3) (2 + v^) (^z^ (2 + yr)' - lot 

2702^/3) -3-5^4/3 



The solution to 



where 



IS 



0(4/3) (2:, t) 



(2 + 77) 
3 ■ 2t 

20 



-$(^,T) 



and 



(29) 



(30) 
(31) 

(32) 



As we will discuss, the universal part of is given by 0(4/3) (z, t). Up to rescaling, the relevant 
behavior of this function is contained in ^{z,T). The higher-order contribution 0(5/3)(;i;, t) 
will be used to compute the two domain function in Section 5.1. 



3.2 Dual theory disk amplitude 

From the dual cubic ([T5|) and the relation g = g/ V2, we can use (^) , (^3|) , (|2^) and (|25|) , 
with 

9c (13 - V7) 

X = — — ^exp(-e2;) , (33) 



to obtain the critical expansion for 0. We find ||T7|| 
2 ( 1 -|- \/7 ) 

= ^ (1 - v^e^ + e^/30(4/3) (^, t) + e^/^0(5/3) (^, t) + 0{e')) , (34) 



with 

and 
where 



<I^W3)i^,t) = ^^^^z,T) (35) 
0(5/3) (2;, t) = , (36) 

(l + v/7) 

Notice that 0(4/3), which expresses the universal behavior in the continuum limit, involves 
the same function $ (given by (53)) which arose in the original Ising theory. 
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3.3 Comparison with Liouville theory 

To compare the universal scaling behavior of the disk amplitude in the Ising and dual theories, 
we recall some of the basic results of Liouville theory. These results were first given in 
23, ER BB, for a review see p^. In Liouville theory the scaling behavior of the partition 



function Z{A) on a sphere of fixed area A is 

Z{A) = A^-^ , (38) 

where for the c = 1/2 model, F = —1/3. On a disk with area A and boundary length I, the 
partition function for conformally invariant boundary conditions scales as 

Z{1, A) = /-3+Q/7^-Q/7e-«VA ^ ^-2/3^-~7/3^-P/A ^ ^39) 

where Q = 7/^6, 7 = ^3/^2. 

The conceptually simplest way to compare the predictions of Liouville theory to the 
matrix model is to directly translate ( ^8|) and (|39D into statements about the asymptotic 
behavior of coefficients in the matrix model generating functions. For example, consider 
expanding the matrix model partition function Z{g) as a function of (7 at c = Cc- 

00 

Zig) = Y: zng"" . (40) 

n=0 

Since n is proportional to the area of the surface, to recover (|38D in the continuum we require 
that Zn — > g~'^n~^'^^^ asymptotically as n — 00. As another example, consider ps, this is 
the partition function for a disk with three boundary edges, which may be thought of as a 
sphere with a single triangle marked. We can expand 

P3ig) = T.^ng'' . (41) 

n 

The coefficients m„ essentially count the total Boltzmann weight of all triangulations with 
n + 1 triangles, with one triangle marked, and should scale as 

nin^n- g^^'n'^^/^ = g-''n~^'^ . (42) 

This scaling behavior can be verified numerically by expansion of the exact solution (^) to 
high order in g. 

An advantage of this way of comparing matrix model results with the predictions of 
Liouville theory is that there is no possible confusion in extracting "universal" parts of the 
matrix model results. The usual approach to this comparison, however, is to take the Laplace 
transform of (^) and (|39D. This approach has the advantage that the calculations can be 
performed analytically; however, we now must be careful in extracting the portion of the 
matrix model result which actually contains information about the asymptotic form of the 
generating function coefficients. A rule of thumb which is usually reliable is that the first 
nonanalytic term in the expansion of an amplitude around e = gives enough information 
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to find the asymptotic form of the coefficients and to check agreement with Liouville theory. 
Let us consider a simple example. It is not hard to verify that when c„ — >■ rfl^~^ as n — > oo, 
with s, t relatively prime integers, we have 

^ c„e-^" = /(e) + r(s/t)e-^/* + 0{e-^'-^^l') , (43) 

n 

where /(e) is analytic in e. The function /(e) is referred to as "nonuniversal" , since it is not 
determined exclusively by the asymptotic behavior of c„, but depends also on the specific 
values of the c„'s for small n. The first nonanalytic term, on the other hand, is determined 
only by the asymptotic behavior of c„, and is thus considered "universal". Thus the rule 
of thumb works in this simple example (although it would fail \i sjt were an integer). In 
Section ^, however, we will find that in some situations a more careful analysis is required. 
We now consider the disk amplitude. Writing as a power series in g 

oo 

Pfc = E S"(lKn , (44) 
n=0 

in order to recover the Liouville prediction (|39|) we require asymptotic behavior 

qk,n - g^-u^'k'/'n-'/'e-''/- . (45) 

(Note that we have an extra power of k since one edge is marked.) To see that this prediction 
is satisfied, we take the discrete Laplace transform 

oo 

J2 e-^'^^-^'"*A;^/3n-'/='e-'='/" = fiez,eh) + 2-^/=^r(4/3)r(-4/3)e^/='$(;2, t/4) + ^(e^/^) , 

n,fc=0 

(46) 

where f{ez,e^t) is an analytic function of e. As in the example (^3]), the universal behavior 



is contained in the first nonanalytic term, which is given in terms of scaling variables {z, t/A) 
by the function $ familiar from the disk amplitudes previously discussed. Thus we see that 
the expansions (^) and ( ^If ) of the disk amplitudes cf) and both have universal behavior 
which agrees with the prediction of Liouville theory for the c = 1/2 model. 



4 Spin/Disorder operators on the sphere 



The scaling dimension of the disorder operator in the dual theory was computed in 
and shown to agree with the KPZ/DDK Liouville results. In this section we compute an 
exact expression for correlation functions of two disorder operators on the sphere in both 
the Ising and the dual theories. We show that the results agree and give the correct scaling 
dimension to the disorder operator. This provides evidence that Kramers- Wannier duality 
is maintained for operators in the bulk after coupling to gravity. 

We compute the correlation function of two disorder operators on the sphere by consid- 
ering a cylinder with a single disorder operator on each boundary; taking the boundaries to 
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have a length which vanishes in the scahng hmit, we obtain the correlation function of two 
disorder operators on the sphere. Notice that, in order to consider boundaries with a single 
disorder operator, it is necessary to introduce a twist in the Z2 bundle over the cylinder 
of which the spin configuration forms a section. (Even so, the total number of disorder 
operators on the entire surface must be even.) There is no way to introduce such a twist 
if the fundamental variables simply represent up and down spins; instead we are forced to 
use variables living between the spin sites, which indicate whether the spins on either side 
are like or unlike. Such variables are natural in the dual theory, but are incompatible with 
the formulation of the original Ising theory as a 2-matrix model. In the calculation below. 



therefore, we use an alternate set of variables introduced in in which boundary data are 
given in terms of spin correlations. 

4.1 Disorder operators in the dual theory 

We begin by looking at the correlation function of two spin operators in the Ising theory or 
equivalently at the correlation function of two disorder operators in the dual theory, 

cDo = ^ (Tr (U-V) ■ Tr {U - V)) = {Tr Y ■ Tr Y) . (47) 

Geometrically, this corresponds to a cylinder for which both boundaries are a single Y 
edge (equivalent in the continuum limit to a sphere with two marked points where disorder 
operators are inserted). 

To compute uq, we find an equation satisfied by the generating function 

00 

cD(x) = ^ (Tr (FX^) ■ Tr y) x'' (48) 

k=0 



in the dual model. Using the formalism of [|T^], an equation for this generating function 



may be derived by considering the result of removing a single Y edge from triangulations 
contributing to (^). We obtain 

xuj{x) = [x(f){x) + 2g{Q{x) - Qq)] , (49) 

so 

(x(l + c) - 2g)Q{x) = x(j){x) - 2^cDo • (50) 

The left hand side of this equation vanishes near x = g = along the curve x = 2g/ (1 + c). 
It follows that _ 

Cooid) = (^) . (51) 
1 + c VI + c/ 



At the critical point, we have 



= . (52) 

1 + c, ^ ^ 
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The expansion for 0(2^/(1 + c)) is just the expansion for (f){x) given in with z = 0. 
Thus, we have 

Oo = (l - ^e^/^t^^ + 0{eH)^ . (53) 

Below we will see that the universal behavior, represented by the non-analytic term of order 
e^/^, agrees with the prediction of Liouville theory. 



4.2 Disorder operators in the Ising theory 

The correlation function of two disorder operators in the Ising theory, or equivalently of two 
spin operators in the dual theory, cannot be expressed in terms of expectation values of the 
matrices U and or X and Y. However, it turns out that we can do the computation using 
an alternative description of the Ising theory, introduced in Ref . []T1| , in which the boundary 
conditions are specified by labeling the vertices between boundary edges as either "stick" or 
"flip" , depending on whether the spins on either side are like or unlike. These variables allow 
us to include configurations with boundary conditions which do not appear in the matrix 
model formulation. 

we derived a generating equation for the stick/flip generating function in 



In Ref. M 



terms of free variables, and this was used to obtain a cubic equation for the homogeneous 
disk amplitude (all sticks on the boundary) that is identical to equation (|^). In a similar 
way, we can write an equation for the generating function 



UJiS 



(54) 



s=0 



where Uk is the weight associated with a surface with two boundaries, one with a single 
vertex marked as a flip, and the other with k stick vertices and one flip vertex. Such 
a configuration involves a single twist in the Z2 bundle over the cylinder, and cannot be 
described as a correlation function of the matrix model operators U and V. The equation 
for uj{s) reads (following the conventions of Ref. ||14|| ) 

\2s^uj{s) + s^u;(s)(0_(s 



SUJiS] 



2 1-c 



+2g{u;{ 



1) 



uJo — suJi) + 2gsuJi + 2s{(j)-{s) + 1] 



(55) 



where 4>-{s) is the homogeneous disk amplitude for a boundary of all stick vertices, and is 
equal to 20(s) — 1 where 0(s) is the homogeneous disk amplitude (|^) of the Ising theory. 
Therefore (EBi) becomes 



u}{s) 

As above, it follows that 



'1 — c)s — s 
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^0(9) 



2s0(s) - gujQ . 
2s(f){s) 2(l-c)s-2^ 



(56) 



(57) 
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where s{g) is defined by the equation 

{l-c)s-s^(l){s)-g = . (58) 

At the critical point, 

(1 - c)Se - ScVc -9c = , (59) 

where Sc = (4 — V7)gc as in (pGl). The expansion for s can be derived by expanding (pH]), 



(60) 



from which it follows that 



4.3 Comparison with Liouville theory 

We can now compare the above results to the predictions of KPZ/DDK based on Liouville 
theory [p4| , p5| , |26| . The scaling dimension of the disorder and spin operators in flat space is 



Aq = 1/16. The two point function of the gravitationally dressed operators on a sphere of 
fixed area A has scaling behavior 

ujl{A) ~ A2-2^Z(A) , (62) 



where A = 1/6 is computed by 
with a satisfying 



A = 1 - - (63) 

7 

Ao = 1 + ^aia - Q) . (64) 

Thus, we expect ojl{A) ~ A~^/^. For the matrix model to reproduce this behavior, we would 
expect that at c = Cc the expansion of ojQijj) 

uoig) = Wnd" (65) 

would have coefficients Wn with asymptotic behavior 

Wn ^ . (66) 

Taking g = gcexp{—te^) as in (p3D , we have the summation 

^o(^) = E ^nd: expi-nte') = a;^ + o;(4/3)e'/'t'/' + 0{eH) , (67) 

n 

where c<j(4/3) is a constant coefficient. This agrees precisely with the universal terms in 
(0) and (0), verifying that the continuum limits of both discrete models agree with the 
predictions of Liouville theory. The agreement between the Ising and dual theories is evidence 
that Kramers- Wannier duality, which exchanges the gravitationally dressed spin and disorder 
operators in the bulk, is a symmetry of the c = 1/2 theory even after coupling to gravity. 
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5 Disorder operators on the disk 



We now proceed to derive analytic expressions for the correlation functions of multiple dis- 
order operators on the disk, both in the regular and dual theories. We begin by using loop 
equation methods to write the exact expressions for the two operator correlation functions. 
We then compute the continuum limit of these expressions, and compare the results to 
Liouville theory. 



5.1 Ising theory 



The generating function for the two domain disk amplitude in the Ising theory, corresponding 
to the correlation function of two disorder operators which separate the domains, is defined 
by 



uU 



vV 



1-uU 1-vV 



(68) 



kA=l 



Using the techniques of |T^, one can derive the following equation for a{u,v) (essentially 
the same equation was derived in [^) 



cj(n, v) 



cgu^ + gv"^ — (1 — c^)uv'^ + v^v"^ u(f){u) + cv(f){v) 
X l^cgu' {(Piv) - 1) + gv' (0(iz) - 1) - (1 

CU(j)y{u) + V(j)^{v) 



—uvpi [v + cu) — guv 



(69) 



I 2 2 
+U V 



(u0(n) + cv(t){v)) (0(m) + 0(w) - 1) 



where 

is related to 0(m) and pi by 

1 



1 oo 
n=l 



u{g-cu)pi-u (2-0(m)) - {g - u){(t){u) - I) 



(70) 



(71) 



Expanding in e and using separate cosmological constants for the two boundary variables u 
and f , 



u 

V 



gc (4 - Vt^ exp{-ezu) 
gc (4 - V7^ exp(-ez„) 



(72) 
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we get 



1 

10 



(9 - 2^7) - (23 + lOv^) e {z,, + z,) 

+9 (3 + V7) e^/=^ (0(4/3) (2:„, t) + 0(4/3) (2;., t) 

+ 9 (3 + Vr) e^/=^ (0(5/3) + 0(5/3) 

+2-8/3 (2 + y7)"/%5/3s(^„,^^,T) + 0{e'] 

where T is given by (|3l|) and the function S is given by 

3T4/3 _ [$(^, T)2 + r)<l>(z, T) + ^y, T) 



z + y 



(73) 



(74) 



We shall see below that J]{z,y,T) plays a similar role for the two-domain amplitude that 
$(z, T) played for the one-domain amplitude, capturing the universal behavior in the con- 
tinuum limit. 



5.2 Dual theory 

In the dual theory, the two-disorder-operator correlation function is given by 



a(x, w) = — (tt ( i Y i Y 



Sk,lU 



k,l=0 



Using the method of ||Tj] it can be shown that a satisfies the equation 



a(x,w) = \ xw {—gx — gw + {1 + c)xw) 

X {^—gx — gw + 2(1 — c)xw — 2gpixw + x^w^0(x)0(w) 
+gw(l){x) — (1 — c)xw(f){x) + x^w(l){xY 
+gx(j){w) — (1 — c)xw(f){w) + xw'^(j){wY 



(75) 



(76) 



As in the previous section, we expand with different cosmological constants for each domain 

g, f 13 - v^) 



X 



w 



6^2 

gc (13 - Vi 

6^2 



exp(^-ez^ 



exp(-ez^) . 



(77) 
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The critical expansion of B{x, w) is then 



a = ^ [l8 (2 + v^) - 18 (I6 + 5V7) e {z, + z^) 

+72 (2 + v^) e^/=^ (0(4/3)(^.,t) +0{4/3)(^«,,t)) (78) 
+3 • 2-^/3 (5 + V7){l + Vrf e^/^S(z., z^, f) + Oie') 

where S(z^, z^, T) is again given by (0). The reappearance of this function is an indication 
of duahty, as we explain below. 

5.3 Comparison with Liouville theory 

To compare these results with Liouville theory, we recall the work of Martinec, Moore, and 
Seiberg in which the scaling dimensions of boundary operators in Liouville theory were 
derived using methods analogous to those of KPZ/DDK. According to their analysis, the 
scaling dimension of a boundary operator with fiat space scaling dimension Am is given 
by 

Afc = 1 - ^ , (79) 
7 

where 

Aw = 1 + «b(2a;, - Q) . (80) 

The flat space scaling dimension of the spin and disorder operators was shown by Cardy to 
be AfjQ = 1/2. Thus, we have 

Afc = 2/3 . (81) 

It follows that the correlation function of two disorder operators on the boundary of a disk 
with area A and boundary length I should be of the form 



R{1', I - I', A) = /-2/3+(l-2A,)^-7/3^ ' ^ ^ ^ (82) 




where /' and /— /' are the boundary lengths separating the operators, and r is an undetermined 
function of the dimensionless ratios. 

Some care is needed to compare (|82| ) with the matrix model results. From the scaling 
dimension of R we see that the universal term in a should scale as e^^^. However, the leading 
nonanalytic term in a is of order e'^^^. In fact, this leading term is not universal, but arises 
from contact terms where the distance between the two disorder operators vanishes in the 
scaling limit. An easy way to see this heuristically is that the universal behavior of a should 
be unchanged if we subtract from it only those terms where the two disorder operators are 
adjacent. This corresponds to taking a{u, v) a{u, v) — 4>v{u) — (pv{v). The expansion of 0„ 
is essentially the same as the expansion of 0, so the effect of this modification is to change 
the coefficients in front of 0(4/3) and 0(5/3) in the expansion ( [73D of cr. 
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This suggests that the only universal terms in a and a are those which are singular in 
both boundary cosmological constants. We now perform a simple calculation which serves 
as a more concrete demonstration of this fact. 

Integrating out A from (^2]), with a vanishing bulk cosmological constant, we have 

S{l\l-l')=l-'"\s[jy^ , (83) 

where s is an undetermined function. The simplest function with the correct scaling and 
symmetry properties is 

S{1', l-l')= /-^/3(;'-4/3 ;/)-4/3) _ ^g4) 

Imagine that the matrix model exactly reproduced this asymptotic behavior; we would then 
have 

~su,i-u-'-\k + l)-'l'{k-"' + l-"'). (85) 
Summing over all boundary lengths, this would give 

oo 

a = Y.e-<'''-+''^\k + lY'"'{k-^'^ + l-^'^) 

k,l=l 

= cr, + cr(i)e(^„ + z,) + cr(4/3)e^/3(zy=^ + zy3) (86) 
+5^(5/3)e^/'(4/' + z'J') + S(^„, z,)e'/' + 0{e') , 

where cr(4/3), o'{5/3) are nonzero constants, and S(z„, z^) is a linear combination of hyper- 
geometric functions. In this case, even though the coefficients / agree with the Liouville 
theory predictions, extra nonuniversal terms are generated which are not simultaneously sin- 
gular in Zu and z^- Thus, in this example, the universal behavior is described by the leading 
term which is nonanalytic in both Zu and z^. 

The hypergeometric function S is not precisely equal to the universal form E, indicating 
that the exact coefficients Sk,i are more complicated than Sk,i- However, from the discussion 



above we see that the universal parts of the expansions (|73D and (|78|) are contained in 
the terms proportional to S, which are singular in both boundary cosmological constants. 
Because these terms take an identical form, we see that the correlation function of two 
boundary disorder operators is identical in the Ising and dual theories. This result clears up 



the puzzle from [13|, and gives strong evidence that Kramers- Wannier duality is preserved 



in the presence of gravity, even on an open string with nontrivial boundary conditions. 



6 Conclusion 

We have calculated the exact two-point functions for spin and disorder operators on the 
sphere and on the boundary of a disk in the c = 1/2 string, both in the Ising and dual Ising 
matrix model formulations. We find that the universal parts of the correlation functions are 
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equal in both models, suggesting that the Kramers- Wannier duality symmetry of the c = 1/2 
conformal field theory is also a symmetry of the full c = 1/2 string. 

In the c = 1/2 CFT, Kramers- Wannier duality is broken by the GSO projection operator 
28| which divides the full Hilbert space into two parts, corresponding to the two dual 



theories. Presumably an analogous procedure is at work in the gravitationally coupled 
theory, although it is unclear how to describe such a mechanism from the point of view of 
the discretized theory. 

Since the duality of a single Ising spin seems to be preserved under coupling to gravity, 
it is interesting to consider the result of coupling two Ising spins to gravity. It is known 
that the CFT describing two Ising spins is a c = 1 orbifold theory with r = 1 [Q. The 
Kramers- Wannier duality of the two Ising spins in this theory gives a Z2 x Z2 symmetry 
to the Hilbert space. It would be interesting to understand the relationship of this duality 
symmetry to the usual c = 1 duality. 

It is straightforward to generalize the methods of this paper to calculate correlation 
functions between arbitrary numbers of disorder and spin operators in either the Ising or 
dual theory. (A calculation along these lines was done in |13|). For a correlation function 



of 2k disorder operators on the boundary of a disk, the Liouville prediction is that the 
correlation function should scale as 7/3 — 2fc/3. Using the methods of |l^, it is also possible 



to calculate these correlation functions in the presence of a boundary magnetic field; this 
work is currently in progress. 
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